clear all; close all; clc;

alpha_f = deg2rad (50); %угол места
 beta_f = deg2rad (0); %азимут
 
%  направление фокусировки
alpha_s = deg2rad (-90:5:90); 
 beta_s = deg2rad (-180:5:180); 
 
 beta = H(alpha_f,beta_f);
 
 F=nan(length (alpha_s),length (beta_s));
 
 for a= 1:length (alpha_s)
     for b= 1:length (beta_s)
         S=1;
         
         y = H(alpha_s(a),beta_s(b))*S; %H=beta
         U = beta'*y;
        F(a,b) = abs(U/S)^2;  % поверхность ДН
     end
 end
 
%  figure;
%  surf(alpha_s, beta_s, F)
 
 [alpha_sm,beta_sm] = meshgrid(alpha_s,beta_s);
 
 
%  surf(alpha_sm, beta_sm, F)
 
[Fx,Fy,Fz] = sph2cart(beta_sm, alpha_sm, F.');  %  перевод из сферических координат
figure;
surf(Fx,Fy,Fz)
 